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Abstract 



In this article we study the dynamics of coupled os- 
cillators. We use mechanical metronomes that are 
placed over a rigid base. The base moves by a motor 
in a one-dimensional direction and the movements of 
the base follow some functions of the phases of the 
metronomes (in other words, it is controlled to move 
according to a provided function) . Because of the mo- 
tor and the feedback, the phases of the metronomes 
affect the movements of the base while on the other 
hand, when the base moves, it affects the phases of 
the metronomes in return. 

For a simple function for the base movement (such 
as 2/ = ^x[^Oi + (1 ~ '")^2] in which y is the velocity 
of the base, is a multiplier, r is a proportion and 
Oi and 62 are phases of the metronomes), we show 
the effects on the dynamics of the oscillators. Then 
we study how this function changes in time when its 
parameters adapt by a feedback. By numerical sim- 
ulations and experimental tests, we show that the 
dynamic of the set of oscillators and the base tends 
to evolve towards a certain region. This region is 
close to a transition in dynamics of the oscillators; 
where more frequencies start to appear in the fre- 
quency spectra of the phases of the metronomes. 



1 Introduction 

Dissipative systems with a nonlinear time-delayed 
feedback or memory can produce chaotic dynamics 
[TJ [5]. The effect of the delay on the dimension of 
these chaotic attractors is shown in [3]. Delay sys- 
tems generically have families of periodic solutions, 
which are reappearing for infinitely many delay times. 
As delay increases, the solution families overlap lead- 
ing to increasing coexistence of multiple stable as well 
as unstable solutions [4j. Anticipating chaotic syn- 
chronization is discussed in [S]. 

Packard showed adaptation to the edge of chaos 
in cellular automata rules with genetic algorithms. 
Some of his results were later disputed in [7]. Co- 
evolution to the edge of chaos is discussed in [S]. 
Edge of chaos has been found to be the optimal set- 
ting for control of a system [S]. A self-adjusting sys- 
tem is a system in which the control of a parameter 
value depends on previous states of the system }10j . 
The authors in [TT] describe adaptation to the edge 
of chaos in logistic map. They believe that adapta- 
tion to the edge of chaos is a generic property of the 
systems with a low-pass filtered feedback. They be- 
lieve that this property is independent of the form of 
the feedback and the system under study. The find- 
ings have also been confirmed experimentally with 
Chua's circuit [IT. In [IS" conserved quantities are 
used for investigating adaptation to the edge of chaos. 
The phrase Edge of Chaos was originally proposed 
by Chris Langton in 1990 in the area of cellular 



automata [13] although some others mentioned the 
same concept around the same time |15j . Guiding an 
adaptive system through chaos is also considered in 
[IB] . Adaptation to the edge of chaos is studied in 
[17j . In that research, one-dimensional chaotic maps 
are used as dynamical systems and feedback of val- 
ues obtained from observations of the system variable 
are provided for the system. The parameter changes 
more slowly than the variable. The authors believe 
that separation of time scales is necessary for the 
method to be self-contained without external control. 
The parameter governs the dynamics of the variable 
every certain number of steps of time series. On the 
other hand, the variables of the system determine 
the dynamics of the parameter. Extinction in a sim- 
ple ecological model is shown in [TB]. It is stated 
that chaotic dynamics does not necessarily lead to 
population extinction. The topic is also discussed in 

[I91E01ET1E2]. 



Purpose of this Research 

We are studying coupled oscillators. We study the 
case in which we interfere with the coupling and for 
this case, we use the concept of forced oscillators and 
use mechanical metronomes placed on a rigid base. 
We study the effects of movements of the base on the 
dynamics of the oscillators. This makes the whole set 
a dynamical system which contains two sub-systems 
(the oscillators) that has many potentials to be ex- 
plored. We will consider the case where the system 
can adjust itself; In other words, we study how the 
system evolves in time when some parameters of the 
system are adjusted by states of the system itself. 



2 Forced Oscillators 

In this section we investigate the effect of forcing the 
oscillators on the dynamics. To read more about 
details of the experiments, please see Appendix [X] 
and to read more about synchronization of oscillators 
when the base is freely moving, please see Appendix 



2.1 Forced Oscillators 

When the metronome base moves, we have a forced 
oscillator for which the governing equation is as below 
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in which x is the horizontal position of the base. We 
write the governing equation as below 

dO 
~dt 



CO 



(2a) 



dio 
Itt 



— —k sin (6*) — e 



- 1 



(2b) 



When we form a feedback loop for the movements of 
the base, we have the following set of equations 
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and h{9,uj) is a function specifying the feedback. 
This function depends on how the base responds to 
the states of the system such as the detected angles 
of the metronomes. For instance, h{9,uj) — 9 means 
that the cart sets its acceleration to the same as the 
phase of the oscillator. 

2.2 Two Oscillators on a Controlled 

Base 

In this case we consider two metronomes placed on a 
moving rigid base. The metronomes oscillate freely 
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while the movements of the base impacts them. The 2.4 Experiments 

equations of motion for this case are as fohows 



LUi = — fci sin {9i) — e 

02 = OJ2 

W2 = -^2 sin (6*2) - e 
i = h{9i,uJi,d2,0J2) 
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in which 61 and ^2 are the angles the metronomes 
make with the vertical and wi and u)2 are their deriva- 
tives (i.e. rotational velocities), x is the horizontal 
position of the base, y — i is the linear velocity of 
the base and z — y is the acceleration of the base (we 
have only included z because it is the only parameter 
affecting the oscillators) . / is the moment of inertia of 
the pendulum, m is the mass of the pendulum, rem. 
is the distance of the pendulums center of mass from 
the pivot point and g is the acceleration of gravity. 

2.3 y = lAre, + {l-r)e2) 

As the servo motor in the experiments receives de- 
sired velocity as the input (and not the acceleration 
or first derivative of acceleration) , we design the func- 
tion for change of z based on a function for y. For 
this case we consider that 



The motor in the experiments receives velocity as the 
desired input and this was the main reason we chose 
y — f {(^1,(^2) format in the model to match the ex- 
perimental limitations as shown in Eq. 
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y = 7.(^^1 + (i-'^)^2) 
z = y 



(9a) 
(9b) 



in which 7^; is an attenuating factor. Fig. [2]shows the 
frequency spectra for experimental results described 
in Eq. "Jp 
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Figure 2: Frequency spectra of 9i — 62 for different 
values of r (ratio) in the experiments. In the experi- 
ments r was fixed and the dynamics was observed for 
the specific value of r. 



y^i.{r9, + {i-r)92) (6) 3 Adaptation 



so 

z = y^-/y{r9, + il-r)92) (7) 

and 

z ^ lz{r9, + {1 - r)92) 
= 7z(rcji + (1 - r)uj2) 

A sample of frequencies for a range of ratios is shown 
in Fig. [1] 



In this section we show how the system evolves when 
we introduce the possibility of adaptation. 

^ Please note that we have ignored the delay between setting 
the desired velocity and the time the move is applied to the 
cart. 

■^To find out more accurate results, we need to run the ex- 
periments for much longer times. On the other hand, there 
are mechanical limitations (such as the limited metronomes 
wind). This is currently a limitation which leaves work for 
future research. 
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Frequency spectrum of 8i for ki = 6.3 and k2 = 6.3 



Frequency spectrum of 02 for ki = 6.3, fe2 = 6.3 




Frequency spectrum of 9i for ki = 6.3, k2 = 18.7 Frequency spectrum of 02 for ki = 6.3, k2 = 18.7 

Figure 1: Frequency spectra for a range of r. In all the experiments 9q — 0.3 and e = 0.64. The vertical axis 
shows the values for the parameter r (ranging from 0.0 to 1.0), the horizontal axis shows the frequencies in 
the signal and the color shows the intensity of the corresponding frequency (blue shows lower intensity and 
red shows higher intensity) 



3.1 Two Oscillators on a Controlled 
Base with Adjusting Parameter 



In this case we also consider two metronomes placed 
on a moving rigid base, but we introduce a param- 
eter for changing the behavior of the system. The 
metronomes oscillate freely while the movements of 
the base impact them. The movements of the base 
are a function of state variables with adjustments 
from a parameter. The equations of motion for this 
case are as shown in Eq. [lO] 
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in which LPF(0i,uji, 62, '^2-, z,p) in Equation lOd 
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low pass filter for modifying the adjusting parameter 
(p). A sample function for LPF is 
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/ = LPF{euUJu02,u:2,z, f)^f- (-0.9/ + O.90i) 

Low pass filter 

(11) 



3.2 Low-Pass Filtered Feedback 

A low-pass filter is made to provide feedback. The 
feedback is used for adjusting the parameter which 
specifies the qualitative behavior of the system. The 
equivalent differential equation is written as 
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So, for our filter, we write the differential equation as 
follows 



dVout 
dt 
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In this filter the cutoff frequency (in hertz) is defined 
as 



fc 



27ri?2C 

or equivalently (in radians per second): 

1 
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The gain in the passband is and the stop-band 
drops off at 6 dB per octave as it is a first-order filter. 
Hence, the differential equations of our system will 
change as follows 



p — a 
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(16) 




p(0) = 0.30 



p(0) = 0.45 



Figure 3: Value of p (the adjustable parameter) in 
time for initial values from 0.00 to 0.45. 




In the filter, the feedback is attenuated based on [10] . 
Also we have considered Ri — R2 = R, the gain of 

the filter is negative (set to be -1 for this case) and Figure 4: Value of p (the adjustable parameter) in 

hence the negated value of p is used in the differential time for initial value = 0.50. 

equations. 
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3.3 Adaptation with the Experimen- 
tal Limitations with Simple Func- 
tion 

As it was mentioned earlier, the motor in charge of 
the cart receives velocity as the desired input. For 
this case we have considered 



y = -/xipOi + {I - p)02) 



(17) 



= -fAOiP - d2P + pOi ~ {p ~ 1)62) (18) 



and 



z = - 02P + 2p(0i - 62) + pOi - pe2 + §2) 

(19) 



3.3.1 / = 7^(-/-(0i-02)) 
Assuming that 
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Figure 5: B\ — B2 when the parameter adapts 
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defines the dynamics of the parameter /, in which 
'jf — a (;^) , now we have 



Figure 6: / when the parameter adapts 

replacing right hand sides of uj yields 

i = 7.(^i(7p7/(-/-(^i-^2))) 
- 02hplfi~f ~ (Oi - O2))) + 27p/(^i - W2) 



P = Ipf 

P = Ipf 

= lplf{-.f -{01-02)) 



(21) 



+ p(-fci sin {61) - e 
+ {\-p){-k2 sin {62) -e 
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(23) 



so, we re- write the relation for z as 

z = lMlplf{-f -{0^-62))) 
-02hplf{-f -{61-62))) (22) 

+ 27p/(wi - UJ2) + pdJi + (1 - P)l^2) 



3.4 Adaptation in Experiments 

In an experiment the parameter was set to be ad- 
justed via a low-pass filtered feedback of 6162 state 
variables. The change of this adjusting parameter in 
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Figure 7: p when the parameter adapts 

time for two different tests are shown in Fig. [H) We 
can find a sense of the meaning of these values for 
the parameter when looking at the frequency spec- 
trum for experiments shown in Fig. [2j 

4 Conclusions and Future 
Work 

In this section we have conclusions of this research 
and will provide some recommendations for future 
work. 

4.1 Conclusions 

In this article we studied the dynamics of two oscil- 
lators that are mechanically coupled through a rigid 
base. The base was controlled with a motor and the 
movements of the base affected the metronome rods. 
We showed that how the functions controlling the 
base, determined the dynamics of the system. The 
function in charge of moving the base had parame- 
ters and we introduced adaptation in the values of 
these parameters. We showed that how the parame- 
ters changed and how was the dynamics of the system 
after change of the parameters. We showed that when 
we introduce adaptation of the adjusting parameter, 
it settles to the direction that the dynamics of the 
system is in a regime between order and chaof[^ 

^In other words, the dynamics of the system goes towards 
a region between lowest and highest number of frequencies in 
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Figure 8: Change of the adjusting parameter in two 
different experimental runs (the values are different 
because of different initial conditions in the experi- 
ments) 

4.2 Future Work 

There are some improvements that are suggested for 
the future work of this research. For analysis and 
numerical simulation, wc considered an ideal model 
of the system in which there are not any delays or 
random terms. An area for future work would be to 
include these terms in the model and investigate the 
system. The similarity to the experimental results 
would show the accuracy of the model. 

Using other alternatives for detecting the phases of 
the oscillators is also a possible improvement for the 
system. This might include using sound detecting 



the frequency spectrum. 
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devices |j and encoders for detecting the angles of the 
metronomes 12 Other devices for sensing can be used 
simultaneously so that the accuracy of the data is 
increased. 

The oscillators are mechanically coupled in this 
system. A possible direction of future research would 
be considering a case in which the coupling involves 
delay. This can be realized when the metronomes are 
positioned on two independent carts and each of the 
carts is capable of moving independently. The same 
can be done for when there are random terms in the 
coupling between oscillators. 

Another possible area for extending the current 
work would be to introduce more oscillators. The 
number of oscillators and their positioning besides 
how they are connected to each other are some pa- 
rameters of the system to be determined. A higher 
number of oscillators will certainly introduce more 
complexity to the system and we expect to see more 
complex (and maybe richer) dynamics. We would 
also recommend using other types of oscillators that 
are not necessarily mechanical (for instance electri- 
cal circuits) . A similar system can be replicated with 
other oscillators and the evolution of the system can 
be studied. 

In this research we showed some cases and consid- 
ered the evolution of the system with the relevant 
parameters and initial conditions. A possible area to 
be investigated is studying the space of parameters 
and the evolution of the system for different regions, 
i.e. what is the role of initial conditions and parame- 
ters of the system (such as dynamics of the cart) on 
the evolution of the system. 
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APPENDIX 

A Experimental Settings 

The experimental settings were used to implement 
the proposed model and validate the numerical re- 
sults for experiments. The description and details of 
these settings are provided in this section. 

A.l Description 

We use mechanical wind-up metronomes as oscilla- 
tors that are mechanically coupled through a mov- 
ing base. The metronomes are Wittner's Super-Mini- 
Taktell (Series 880) and are claimed to be the world's 
smallest pendulum metronomes. Another part of the 
system is in charge of reading the angle values of 
metronome rods in time. A schematic of the system is 
shown in Fig. [9] and some images of the experimental 
settings are shown in Fig. 10 




of the metronomes and hence might have affected 
the dynamics, we preferred to choose an approach 
for reading the values of angles that does not affect 
the dynamics of the system. 

We selected a vision and an image processing ap- 
proach so that we track the rod; for instance, the 
conservation of mechanical energy using video anal- 
ysis has been investigated in [53]. We used simple 
colored stickers so that the vision system can use the 
color for tracking the position. The stickers are cir- 
cles for which the radius is 9.4 mm and the choice of 
colors are red, green, yellow and blue. Regarding the 
position of the stickers, the center of each sticker was 
marked and a tiny hole was made in the center. This 
hole was used to align the position of the sticker on 
the center of the rod. 

A. 3 Image Processing 

The main image processing library that is used in 
this system is OpenCV |^ A simple USB webcam 
is used to capture the image and the capture frame 
rate is around 20 to 40 frames per second We do 
not use a pre-recorded stream for processing so that 
we can have real-time control on the system. In or- 
der to track the objects in real-time, an algorithm 
is in charge of looking for proper pixels in a cer- 
tain window around the object that is being tracked. 
With enhanced settings, we used more than one color 
sticker for each metronome rod so that we increase 
the accuracy of angle detection (color image process- 
ing has been used in some applications such as [25]). 
A sample of the perceived image is shown in Fig. Ill] 



Figure 9: The schematic of the system 



A. 2 Settings 

One of the main factors we had in mind was that 
to build a robust and reliable setting that is easily 
available and possible to use. We needed to read the 
values of angles at different times and one of the first 
candidates was using encoders for reading the values. 
As the encoders might have affected the mechanics 



A. 3.1 Finding the Angle 

When using more than two stickers on the rod, we 
need to find the least squares regression line. We 



^OpenCV is a computer vision library originally devel- 
oped by Intel. It is free for use under the open source 
BSD license. The library is cross-platform and it focuses 
mainly on real-time image processing. It is available at 
http:/ /opencv. willowgarage.com/ 

^This is a rough estimate and is extracted based on exper- 
iments in normal conditions. It is changed according to the 
settings of the program, hardware specifications and the cam- 
era used for image processing. 
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Stickers used for image pro- Two metronomes on a freely Servo motor and the controller A general view of the whole 
cessing moving cart for moving the cart settings including the camera, 

the track, the carts, the motor 
and the pc 

Figure 10: Experimental Settings 




Figure 11: Two metronomes on two connected carts. 
The crosses show the center of mass of the detected 
stickers. Also there are rectangles around each of the 
stickers which show the area to be searched at each 
step. 



use the following formula for finding the angle of the 
metronome rod 



/ n n n \ 

E - E E 

i—l i—1 i — 1 



tan 



n E 

1=1 



E 

1=1 



(A-l) 



in which n is the number of stickers and Xi and yi are 
the X and y of the ith sticker in pixels. 



A. 3. 2 Validating the Results 

The image processing module might occasionally lose 
track of stickers. To prevent the consequences of this 
problem and to validate the results, correlation of 
the positions of the detected stickers are calculated 
as below 



xy 



(A-2) 

in which n is the number of stickers and Xi and yi 
are the x and y of the ith sticker in pixels. If corre- 
lation is more than a certain threshold, the sample 
is accepted, used for updating the values and stored. 
If the correlation is smaller than the threshold, then 
there might be one possibility: The positions of the 
detected stickers are too close to a vertical line (which 
reduces the correlation). In that case, the detected 
points are virtually rotated 7r/4 radians and the cor- 
relation is calculated again. If the correlation is more 
than the threshold, it means that the detected points 
have been sufficiently close to a line and the line has 
been vertical on the first round. If correlation is less 
than the threshold even in the second round, then it 
shows error in detecting the stickers and the sample 
is discarded. 
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A. 4 Model of Oscillators 

The equation that governs the motion of a single 
metronome is described as [231 IMl 



(P9 mrc.m.g 




sin {9)+e 



- 1 



de 

di 



= (A-3) 



which is rewritten as 



Lu = — fc sin (9) — e 



- 1 



(A-4a) 
(A-4b) 



in which 9 is the phase of the metronome (angle made 
with the vertical line), / is the moment of inertia of 
the pendulum, m is the mass of the pendulum, rem. 
is the distance of the pendulum's center of mass from 
the pivot point, e specifies the effect of escapement 
and damping, is the van der Pol term and g is the 
acceleration of the gravity, k — ""'y determines 
the frequency of oscillation. The damping function 
D{9) is defined as ^ 



m = 



- 1 



(A-5) 



and with small e, this term produces small oscillations 
with an amplitude of approximately 29q in an isolated 
oscillator 23; . 

To find out the parameters based on experimental 
data, we set a metronome to different frequencies and 
let it oscillate for about 1000 seconds. The data was 
collected as pairs of time and angle and used for sys- 
tem identification. Using the same initial conditions, 
a numerical simulation was run with some initial pa- 
rameters. The error was defined as 



Error 



N 

i=i 

N 



(A-6) 



in which 9{t) shows the value of angle in the simu- 
lation at time t, S[i, 1] is the time of capturing the 
iih sample, S[i,2] is the angle value of the ith sam- 
ple and N is the total number of experiment samples. 



The identification method was to minimize the error 
between the experimental data and numerical sim- 
ulation of the differential equations. At each step, 
two of the parameters were fixed and one parameter 
was modified in the direction of minimizing the er- 
ror. When the error was minimized, the parameter 
was fixed and another one was modified. The pro- 
cess was repeated for several iterations until all the 
parameters settled in a range that was sufficiently 
small. For this case we identified the parameters k, 
e and 9q. Fig. [T2| shows the error value for each of 
the parameters when the other two parameters are 
fixecFl Fig. 13 shows the experiment samples and 



the proposed model and Fig. [Mjshows the error when 
changing the parameters in iterations. 



Figure 14: Identification error when one of the pa- 
rameters is changed during iterations 



"Naturally the best value of the parameter on the horizontal 
axis is the value for which the error on the vertical axis is 



minimum. 
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Figure 12: Identification Errors for three parameters 



9 (radians) 




Time (seconds) 



Figure 13: Identification of the model: The continuous graph shows the numerical simulation of the proposed 
model and the dots show data extracted from the experiments 
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B Synchronization 

Synchronization is a ubiquitous phenomenon every- 
where and it plays an important role in many phe- 
nomena. We start by reviewing some literature about 
oscillation, coupled oscillators, synchronization and 
some related areas. 

B.l Oscillation and Coupled Oscilla- 
tors 

The concept of self-oscillations was first proposed by 
Andronov, Khaikin and Vitt in 1937 ^ (for the En- 
glish version see [28] )• A self-oscillating system is 
considered as "an apparatus which produces a peri- 
odic process at the expense of a non-periodic source 
of energy." Appleton ^ and van der Pol [501 131] 
set experiments with electric circuits while they were 
studying the reception of radio signals with electric 
circuits with triodes. In a relaxed self-sustained os- 
cillator, although the motion of a point in the phase 
plane might be non-uniform, but the growth of the 
phase in time is still uniform |32j . 

The interaction between the organ pipes was stud- 
ied by Rayleigh [33! • His experiment involved two or- 
gan pipes which had close peach and were placed close 
to each other. The pipes together sounded in per- 
fect unison. Another investigation is done in [53] in 
which one pipe was substituted by an electric speaker. 
The authors observed that even minute driving sig- 
nals forced the pipe to synchronization. 

B.2 Synchronization 

The word synchronous originates from the Greek 
words chronos and syn which means "sharing the 
same time" |32j in which synchronization is consid- 
ered as a complex dynamical system rather than a 
state. Synchronization is discussed in j35| . It can 
result from an interaction between systems [53] or 
sub-systems [5711551155]. 

Synchronization of two coupled escapement-driven 
pendulum clocks was investigated in [40]. Blekhman 
[4T] discusses observations of the Dutch researcher 



Christian Huygens [j and presents the results of a 
laboratory reproduction and a theoretical analysis 
of oscillators coupled through a common supporting 
frame. At first, Huygens suspected the "sympathy" 
between the clocks was due to induced air currents, 
but eventually concluded that the cause was the "im- 
perceptible movements" of the common supporting 
structure [42]. In the reproduction of Huygens ex- 
periment, the anti-phase state was the only type of 
synchronization that was observed [32]. 

The first observations of synchronization in elec- 
tronic tube generators were done by Eccles [351 133] ■ 
He considered the problem of creating a precision 
clock and the transmission of naval signals. In [3S] 
the authors experimentally studied rhythmic hand 
clapping. They consider a conflict between average 
noise intensity and synchronization. They discuss a 
mechanism of hand clapping period doubling by indi- 
viduals that helps the group achieve synchronization. 
Their results offer a novel route to synchronization, 
not observed in physics or biological systems by the 
time of publication. 

Synchronization only happens in self-sustained sys- 
tems [46 . The authors describe self-sustained oscilla- 
tors mathematically as an autonomous (i.e. without 
explicit time dependence) nonlinear dynamical sys- 
tems. Phase of an oscillator is considered neutrally 
stable. This means that a slight perturbation can 
change the phase while the amplitude is stable and 
is not affected by external perturbations. This prop- 
erty provides the oscillator the ability to synchronize 

m- 

Synchronization of two metronomes is investigated 
in [47] . The authors address the problem of analytical 
study of in-phase synchronization for the model of 
two metronomes an the common support proposed 
in [55]. The authors proved the existence of an in- 
phase regime and they proved that for when the angle 
difference of metronomes was zero (0i — 02 =0), the 
sum of angles {4>i + ^2) has a periodic regime. 

The synchronization of an array of clocks hang- 
ing from an elastically fixed horizontal beam is stud- 
ied in [3H]- The beam is considered as a rigid body 



^Christian Huygens is probably the first scientist who ob- 
served and described the synchronization phenomenon. 
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connected to a spring and a damper. Different 
types of synchronization were observed: Symmetrical 
Synchronization, Complete Synchronization and De- 
synchronous Behavior. Synchronization of two and 
more metronomes is considered in . Synchroniza- 
tion of coupled mechanical metronomes is also stud- 
ied in [22]. The authors study synchronization by 
means of numerical simulations showing the onset of 
synchronization for two, three and 100 globally cou- 
pled metronomes. 

Crowd synchrony on London Millennium Bridge 
P5] is investigated in [5D]. Sufficient conditions for 
controlled synchronization of non-linear systems is 
provided in |51j . An attempt to provide a general 
formalism for synchronization in dynamical systems 
is shown in |52) . Frequency and coordinate synchro- 
nizations are considered in that article. It is men- 
tioned that synchronization as a phenomenon should 
be considered in context and depends on the view. A 
system showing synchronization viewed from a view 
point, might not seem having synchronization from 
another point of view. 

A sample of synchronization is observed in various 
species. For instance, banded mongoose groups show 
high degree of birth synchrony to avoid the negative 
effects of competition with other females [S3]. Syn- 
chronization of time-delayed systems is discussed in 
[54] . A brief introduction to the theory of synchro- 
nization of self-sustained oscillators is presented in 



B.3 Model 

A simple model of coupled oscillators is given by 



9i =uji+Ki sin {62-01) 
h =UJ2+K2 sin (01 - 62) 



(A-7) 



in which 61 and 82 are the phases of the oscillators 
and uii and W2 are the natural frequencies. The two 
parameters Ki and K2 determine the amount of de- 
pendency between the two oscillators (and show that 
how much they are affected from the other one). For 
an uncoupled system, we have Ki = K2 = and 
hence the equations are reduced to 



i^2 



(A- 



When the oscillators are coupled, there is a phase 
difference 6 between the two oscillators 



(A-9) 



So, we have 



= 02-01 

= uji - UJ2 — {Ki + K2) sin ( 



(A-10) 

in which loi — LJ2 is also called Frequency detuning 



Synchronization is investigated in other fields such 
as Small- World networks [SHI [ST] |SS] , weighted com- 
plex networks |591 1601 161) and dynamical networks 
[52] . Phase synchronization of weakly coupled self- 
sustained chaotic oscillators is investigated in [55] . 
The exact mechanisms of generation of epileptic 
seizures in human brains is still uncertain. Neverthe- 
less, it is widely accepted that an abnormal synchro- 
nization of firing neurons causes epileptic seizures. 
To investigate this problem, phase synchronization 
between different regions of the brain is measured 

IMESl. 



Different settings and sets of parameters were used 
to study synchronization. Some of them and the re- 
sults are described in this section. 



B.4 Oscillators on a Freely Moving 
Base 

When the oscillators are placed on a freely moving 
base(a sample is shown in Fig, 
possibility of synchronization, 
settings is shown in Fig. [l^ 



10 1, they have the 



A schematic of the 
The equations of motion 
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Figure 15: The schematic of the system for two 
metronomes placed on a freely moving base. The 
metronomes move the base and the movements of 
the base affect the metronomes. The camera and the 
PC are only for observation. 



in this case are as follows 



(1 + A) sin (6*1)+/^ 
-/3cos (^'i)^(sin( 
(1- A)sin(02)+Ai 



-/3cos (6*2 )^ (sin ( 



^ -1 



d7 



sm ( 

2 

^1 -1 



0) = o 

dfh 
dT 



l) + sin( 



h)) = 

(A-11) 



in which t — ujt is a dimensionless time variable, 
^2 _ rnr^^^rn.g ^j^g squarc of the average angu- 
lar frequency of the uncoupled, A k, "^""^ is the 
relative frequency difference between the oscillators, 
small amplitude oscillator without damping or driv- 
ing, 00 is the van der Pol term and /3 is the coupling 
parameter 



(A-12) 



M + 2m 



in which M is the mass of the base, / is the moment 
of inertia of the pendulum, ni is the mass of the pen- 



dulum (masses of the pendulums are considered to 
be the same) and rem. is the distance of the pendu- 
lum's center of mass from the pivot point. For low 
damping of the coupling medium we may observe in- 
phase synchronization and when the damping of the 
coupling medium is high, we can observe anti-phase 
synchronization. 

B.5 Synchronization on a Freely Mov- 
ing Base 

The metronomes were set to various frequencies and 
then placed on the freely moving base. The phase 
differences were observed for each frequency. The 
experimental results are shown in Fig. |16[ 
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Figure 16: Phase difference of two metronomes 
placed on a freely moving base in the experiments 
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